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Abstract 

Comparability graphs are the undirected graphs whose edges can be directed so that the 
resulting directed graph is transitive. They are related to posets and have applications in 
scheduling theory. This paper considers the problem of finding a transitive orientation of a 
comparability graph, a requirement for many of its applications. A linear-time algorithm is 
presented based on an elegant partition refinement scheme developed elsewhere for the problem. 
The algorithm is intended as a simpler and more practical alternative to the existing linear¬ 
time solution, which is commonly understood to be difficult and mainly of theoretical value. 
It accomplishes this by using Lexicographic Breadth-First Search to achieve the same effect as 
produced by modular decomposition in the earlier linear-time algorithm. 


1 Introduction 

Comparability graphs are a well-studied class of graphs related to posets. Every poset can be 
represented by a transitive, directed, acyclic graph; comparability graphs are formed by ignoring 
the directions on these edges. Equivalently, comparability graphs are those undirected graphs whose 
edges can be directed so that the resnlting directed graph is transitive: whenever there is an edge 
from a to 6 and one from b to c, there is also an edge from a to c. Given a comparability graph, the 
transitive orientation problem asks you to find such a transitive orientation of its edges. An example 
of a comparability graph is provided in fignre while the same graph with its edges transitively 
oriented appears in figure Graphs whose complement is a comparability graph are called co¬ 
comparability graphs. A related problem is to compute a transitive orientation of the complement 
of a co-comparability graph. Permutation graphs are those that are both comparability graphs and 
co-comparability graphs; these correspond to partial orders of dimension-2. (See [T] for a summary 
of these graph classes.) 

Many interesting combinatorial problems can be efficiently solved for comparability graphs and 
co-comparability graphs (and thus permutation graphs) once a transitive orientation is known jl4j . 
This is the source of their many applications in scheduling theory [HI 13]. Earlier solntions to the 
transitive orientation problem ran in time O(n^) [21] and 0{6m) |9l[8l[T9|, where 5 is the maximum 
degree of a vertex in the graph. The approach developed in the O(n^) algorithm of [2T] was later 
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extended to 0(mlog n) [16] and 0(n + m) |15j algorithms. While the former presents an elegant 
solution that is straightforward in its understanding and implementation, linear-time in the latter 
was achieved at the expense of such simplicity. Although the latter’s algorithm’s difficulty has been 
acknowledged by the authors in |16[ [22] , it has so far remained the only linear-time solution to the 
transitive orientation problem. Part of the acknowledged difficulty is its reliance on the modular 
decomposition algorithm also developed in m- 

Modular decomposition has long been known to simplify the computation of transitive orienta¬ 
tions |7] in addition to numerous other problems [18]. The algorithm in [T3] and another indepen¬ 
dently developed in [2] were the first linear-time modular decomposition algorithms. Unfortunately, 
both these algorithms were complex enough to be deemed mainly theoretical contributions [2] . But 
the practical importance of modular decomposition led many to attempt simpler, linear-time alter¬ 
natives |12L [5] (note that while [10] claimed linear-time, the authors have since noted an error in 
the paper ffH), something finally achieved in [23] by unifying earlier approaches. Similar attempts 
at a simpler, linear-time alternative to the transitive orientation algorithm of |15j have so far not 
appeared. 

The elegant 0(m log n) transitive orientation algorithm of |16j is based on partition refinement 
techniques (see m for an overview). Some hope was expressed there that a simpler linear-time 
implementation of these techniques may be possible. This paper posits an algorithm that proposes 
to answer that question in the affirmative. Where previous algorithms used modular decomposi¬ 
tion, this paper achieves the same effect with Lexicographic Breadth-First Search (LBFS) (see [3] 
for a survey). The key point is that LBFS is itself easily implemented in linear-time using partition 
refinement |12) . The result is a simpler, linear-time alternative to solving the transitive orientation 
problem. The paper also suggests how its approach with LBFS can be adapted to compute a 
transitive orientation of the complement of a co-comparability graph in linear-time, which neces¬ 
sarily avoids explicitly computing that complement. All of this can be implemented using partition 
refinement techniques that amount to basic list and tree traversals. 

2 Background and Overview 

2.1 Preliminaries 

Only simple graphs are considered in this paper. All graphs should be assumed to be undirected 
unless otherwise specified. The set of vertices of a graph G will be denoted V(G) and its set of 
edges E{G). Throughout this paper n will be used to refer to |U(G')| and m for \E\. The graph 
induced by the set of vertices S C V{G) will be denoted G'[S']. The set of neighbours of a vertex x 
will be denoted N{x), while N[x] = V(x)U{x}. Depending on the context, xy will either represent 
the undirected edge {x,y} or the directed edge {x,y). The connected components of a graph will 
simply be referred to as its components] the connected components of the complement of a graph 
will be called its co-components. 

A vertex x is said to be universal to a set of vertices S' if x is adjacent to every vertex in S, and 
X is isolated from S if it is adjacent to no vertex in S. A set S is said to be universal to another 
set S' if every vertex in S is adjacent to every vertex in S' (i.e. there is a join between the two 
sets). On the other hand, x splits S if it is adjacent to at least one but not every vertex in S, and 
a set S splits another set S' if S n S' / 0 and S' — S ^ Meanwhile, two sets S and S' are said 
to overlap if S — S', S n S', and S' — S are all non-empty. A set S C V{G) is trivial if S = V{G) 
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Figure 1: A comparability graph - vertices are grouped in boxes, and edges from one vertex to the 
perimeter of a box means that vertex is universal to the vertices inside the box. 


or |5| = 1. 

The vertices of a rooted tree will be referred to as nodes, with the non-leaf vertices specially 
designated as interior nodes. An ancestor of a node x is a vertex other than x on the path from 
that node to the root. A node x is a descendant of the node y ^ x \i y appears on the path from 
X to the root. 

2.2 Comparability Graphs 

A directed graph is transitive if whenever there are edges ab and be there is also the edge ac. Vertices 
whose incident edges are all directed the same way are either sources (directed outward) or sinks 
(directed inward). An undirected graph whose edges can be transitively oriented is a comparability 
graph. A linear extension of a comparability graph G is an ordering of V{G) that induces the 
following transitive orientation: the edge ab is directed toward b if and only if a appears before b 
in the ordering. An example of a comparability graph appears in figure while an example of a 
transitive orientation of its edges appears in figureObserve that x, z, q, w, r, v, y, u, a, 6 is a linear 
extension for this graph. 

2.3 Modular Decomposition 

A module is a set of vertices M such that no vertex x ^ M splits it. A module M is strong if 
it overlaps no other module. A graph is prime if it contains no non-trivial module. Note that 
the graph in figure is prime. The modular decomposition tree for the graph G is the rooted tree 
defined as follows. 

• The set of leaves corresponds to the set of vertices V(G). 

• The set of interior nodes corresponds to the set of strong modules. 
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Figure 2: The comparability graph of figure with its edges transitively oriented. 


• There is an edge between two nodes if one is the smallest strong module containing the other. 

Every interior node in the modular decomposition tree defines a quotient graph: the graph 
induced by the set of vertices formed by selecting a single leaf descending from each of its children. 
It can be shown that the quotient graph is invariant under different choices of leaves in that the 
resulting graphs are all isomorphic. Every quotient graph is either a clique, independent set, or is 
prime (see, for example, a)- 

Modules play an important role in the transitive orientation of comparability graphs. Given 
a module M in a comparability graph G, and a vertex x G M, it suffices to separately compute 
a transitive orientation for G[M] and G[{V{G) — M) U {x}]: merely replace x in the latter with 
M. This can be applied to the modular decomposition tree’s strong modules and quotient graphs. 
Since cliques and independent sets can trivially be oriented, the problem therefore reduces to the 
prime case. This is made somewhat easier by the fact that prime comparability graphs have a 
unique transitive orientation (up to reversal of the directions on all edges) [7] . 

2.4 Partition Refinement 

A partition of the set of vertices V(G) is a collection of disjoint subsets P = {Pi,..., Pk} such that 
V{G) = Pi U • • • U Pfc. Each Pi is called a partition class. A refinement of the partition V is one 
in which every partition class is a (not necessarily proper) subset of a partition class in V. The 
partition of S' C F(G') induced by P, denoted P[S], is obtained by taking the intersection of S 
with every partition class in V (ignoring any empty intersections). Partition refinement frequently 
operates on ordered partitions: those having an ordering imposed on their partition classes, written 
P = Pi,..., Pfc. 

Refinements in this paper will be computed according to the pivot operation defined by algo¬ 
rithm The input vertex p is referred to as a pivot. Algorithm subdivides partition classes 
according to N{p) as follows: 
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Remark 2.1. Suppose that A,B are partition classes (in order) that replace a partition class P 
during algorithm^ Then A,B^%,p^P, and if p is in an earlier partition class than P, then 
Ar\N(j)) = 0 and B C N{p), otherwise A C N{p) and Br\N{p) = ^. Moreover, if P is not replaced 
during algorithm^ then either p £ P, P N{p) = %, or P C N{p). 


Algorithm 1: Pivot{V,p) jUj 

Input: An ordered partition V = Pi,... ,Pk and a distinguished vertex p ^ Pi. 

Output: A refinement of V. 

foreach Pj ^V,j i such that N{p) splits Pj do 

A^Pj-Nip); 

B Pj - A; 

if i < j then replace Pj in V with A, B in order else replace Pj in V with B, A in order 

end foreach 

return V] 


Observe that in algorithmic it may happen that one of A or R is the empty set. By convention, 
throughout this paper, we take the approach that such an empty set is not inserted into the partition 
as a partition class. Later in this paper, we will have need to generalize the above pivot operation. 
There too we adopt this convention. In general, whenever partitions and refinement are concerned 
in this paper, the empty set is not permitted. 

The pivot operation above can be used to compute a source vertex in a prime comparability 
graph according to algorithm [13] . If the partition class P at the end of algorithm [^ did not 
contain a single vertex, then it would be a non-trivial module, by remark 2.1, contradicting the 


input graph being prime. The fact that this vertex must be a source (or equivalently, a sink) vertex 
is due to the following easily verified invariant: throughout the algorithm, all edges between any 
vertex in P and it neighbours outside P must be directed the same way. Applying algorithm[Cto the 
graph in figure [T] where the sequence of pivots chosen during its execution is a, z, b, u, y, q, w, r, v, x, 
results in the vertex x correctly being returned as a source vertex. 


Algorithm 2: Source{G) [Hj 

Input: A prime comparability graph G. 

Output: A source vertex in the transitive orientation of G. 

pick some vertex x G V{G); 

let V be the ordered partition {x}, V{G) — {x}; 

let P be "P’s last partition class (throughout what follows); 

while there is a vertex z ^ P that has not been pivot do Pivot{V, z) 

return P; 


An ordered partition V oi V{G) is consistent with a linear extension if there is a way of per¬ 
muting the vertices in each of its partition classes that results in a linear extension. Observe that 
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if X is a source vertex, then {x},V{G) — {x} is consistent with a linear extension. Also observe 
that if the input to algorithmis consistent with a linear extension, then so too will be its output. 
Finally, observe that any non-trivial partition of a prime graph contains a partition class P that is 
split by some vertex y ^ P. 

The above observations suggest the following simple algorithm for computing a linear extension 
of a prime comparability graph G: hrst compute a source vertex x using algorithmj^ then repeatedly 
apply algorithmic starting with the ordered partition {x},V{G) — {x}, until all partition classes 
are singleton sets. This approach was first articulated in m and later generalized to the non-prime 
case in |12) . The latter uses the following data-structure to represent ordered partitions: 

• The ordered list of partition classes is represented as a doubly-linked list. 

• The elements within each partition class are represented as a doubly-linked list. 

• Each partition class maintains two indices indicating its range: the number of elements in 
partition classes before it, plus the number of its own elements. 

• Each partition class maintains a pointer to the first element in the doubly-linked list repre¬ 
senting its elements. 

• Each element within a partition class maintains a parent pointer to its containing partition 
class. 

The corresponding implementation of algorithm [C proceeds as follows: 

• Traverse the adjacency list of the pivot to determine the partition classes different than its 
own that contain at least one of its neighbours. 

• Use the range indices to determine the relative position of the partition class containing the 
pivot and those encountered above. 

• For each partition class encountered above, create an empty partition class immediately before 
those that appear earlier than the pivot, and a new partition class immediately after those 
that appear after the pivot. 

• Traverse the adjacency list of the pivot once more, and move each neighbour to the adjacent 
empty partition class created above, updating the range of the affected partition classes to 
reflect this change. 

Algorithmclearly runs in time 0{\N{p)\) using the above implementation. The efficiency of 
the approach described above for computing linear extensions therefore depends on the number of 
times algorithm is invoked. To limit the number of times a vertex is a pivot, |14] developed an 
ingenious rule: only reuse a vertex as pivot when its containing partition class is at most half the 
size it was the last time it was pivot. The result is an 0(m logn) transitive orientation algorithm 
for prime graphs. The same rule was applied in [12] for the non-prime case, the result also being 
an 0(m log n) algorithm. 
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2.5 Lexicographic Breadth-First Search 

In Breadth-First Search, the order in which a vertex is explored depends only on its earliest explored 
neighbour. Lexicographic Breadth-First Search (LBFS) goes further by using all previously explored 
neighbours to determine when a vertex is explored. This is accomplished through the lexicographic 
labelling scheme specified in algorithm throughout the algorithm vertices carry lexicographic 
labels that have been assigned to them by their previously explored neighbours - vertices with 
earlier explored neighbours will have lexicographically lager labels - and on each iteration of the 
algorithm, a vertex with lexicographically largest label is chosen. It is important to note that on 
any iteration of the algorithm, many vertices may share the lexicographically largest label. 

To understand the operation of algorithm consider its operation on the graph in figure 
One possible output of the algorithm is x, b, y, u, z, q, w, r, v, a. To see why, notice that initially, 
all vertices share the same label, and therefore any one of them can be chosen, say x. From 
there, x labels its only neighbour, b, and all other vertices retain their empty label. Therefore b 
has lexicographically largest label and must be chosen. From there, b labels its neighbours, while 
others retain their empty label, and so on. 


Algorithm 3: LBFS{G,x) jJO] 

Input: A graph G. 

Output: An ordering a of V{G). 

initialize a as the empty list; 

initialize each vertex in V(G) by assigning it an empty label; 
while |cr| 7 ^ |F(G)| do 

select a vertex x ^ a having lexicographically largest label and append it to ct; 
foreach vertex y G N{x) — a do append |F(G)| — |cr| -|- 1 to y’s label 

end while 

return cr; 


As we saw above, the initial vertex of any lexicographic breadth-first search is selected arbitrar¬ 
ily. A variant of algorithm allows this initial vertex to be specified as an input to the algorithm. 
This is accomplished by assigning the input vertex some large initial label, say 00 prior to any iter¬ 
ation of the while loop. In proving results and defining things about LBFS we will use algorithm 
as it is specified above. However, it will be convenient throughout this paper to sometimes refer to 
the variant taking an input vertex, with the intention being clear from context. 

Although algorithm produces an ordering of vertices, more important for this paper will be 
the labels it assigns in doing so. These are important in so far as they define slices: any set of 
vertices sharing the lexicographically largest label immediately before an iteration of the while-loop 
in algorithm]^ Hence, there are |F(G)| = n different slices, one for every iteration of the while-loop, 
and the set of slices can be ordered according to the iterations on which the slices are defined. The 
vertex selected from the slice and added to a on each iteration is the slice’s initial vertex. Notice 
that the set V{G) is always the first slice. Building on our example earlier, one possible ordered 
set of slices and initial vertices for algorithm on the input graph in figure is the following: 

• Si = {x, b, y, u, V, r, q, z, tc}, initial vertex x; 


7 






• ^2 = {6}, initial vertex b] 

m Ss = {y, u, V, r, q, z, w}, initial vertex y; 

• Si = {u, V, r, q, z, r/;}, initial vertex u; 

• 85 = {z, q, te}, initial vertex z; 

• Sq = {g, t(;}, initial vertex g; 

• St = initial vertex w; 

• Sg = {r}, initial vertex r; 

• Sg = {?;}, initial vertex v] 

• •S'lo = {a}) initial vertex a. 

It is not difficult to see the hierarchical relationship between slices in the example above. This 
is formalized in the following remark: 

Remark 2.2. If S and S' are two distinct slices, then either S C S', S' C S, or S (1 S' = 0. 

This leads to the idea of maximal subslices: S' is a maximal subslice of 5 if 5' C 5 and there 
is no other slice S" such that S' C S" C S. By convention, the initial vertex of a slice is also 
considered one of its maximal subslices. Further, throughout this paper, we will generally use 
Si,... ,Sn to denote the set of all slices, and x, Si,..., Sk to denote a set of maximal subslices, 
where x is the initial vertex. 

Observe that the set of maximal subslices is clearly a partition of that slice. And of course, the 
ordering of all slices defines an ordering of any set of maximal subslices (with the convention that 
the initial vertex always appears first). 

An edge is active for a slice S if its endpoints reside in different maximal sublices of S. The 
notion of active edges is adapted from [5] , where it was introduced for the hierarchical set of strong 
modules that arise in modular decomposition. Their reason for doing so relates to the following 
observation, presented here in terms of slices. 

Remark 2.3. For any edge there is a single slice for which it is active. 

The following are properties of active edges specific to the LBFS context, each a direct conse¬ 
quence of the labelling scheme employed by algorithm 

Remark 2.4. Let S be a slice and x. Si,..., Sk its maximal subslices in order. Then: 

1. Either N{x) D S = Si or N{x)n S = 0; 

2. Every vertex y G Si is either universal to or isolated from Sj,i < j; 

3. For every Si and Sj, 1 < i < j, there is a vertex y G Si, i < i, such that y is universal to Si 
and isolated from Sj. 


The last two remarks form the basis for the transitive orientation algorithm presented in this pa¬ 
per. They imply many other properties of active edges that will be presented as needed throughout 
the paper. One such property is the following: 


Lemma 2.5. Consider an LBFS ordering of a comparability graph G. Let S he one of the slices 
and X, Si, ..., Sk its maximal sub slices in order. Assume that N{x)r\S / 0. Then in any transitive 
orientation, the edges between all vertices in some Si,i > 1, and any vertex z ^ Si are all directed 
the same way. 


Proof. Assume that z ^ Sj. If j < i, then the result follows by condition 3 of remark 2.4 So assume 
that j > i. Then by the same remark, there is a vertex q £ Si,i < i that is universal to Si but 
isolated from z. As with z, all edges from q to Si must be directed the same way. Therefore all 
edges between Si and z must also be directed the same way. □ 


2.6 Overview 

This paper develops an algorithm to compute a linear extension of a prime comparability graph. 
It uses the same basic approach as was described earlier in the context of partition refinement. 
The difference is that instead of limiting the number of times a vertex is pivot, this paper limits 
what edges incident to a pivot are used each time. The set of edges incident to each vertex will 
be partitioned, and every time a vertex is pivot, a unique subset of these incident edges will be 
processed. In the end, each edge will only be processed once, allowing for linear-time transitive 
orientation, 0{n + m). 

Partition refinement will remain the primary mechanism for computing the linear extension. 
Algorithmic will need to be generalized so that only a subset of vertices in a pivot’s neighbourhood 
is processed. The partition of the set of edges incident to each vertex will be defined in terms 
of the slices for which each edge is active. This will necessitate computing an LBFS ordering, 
the corresponding set of slices, and the active edges they define. The co-components of the graph 
induced by each slice will also need to be computed and treated as a kind of pseudo-slice for reasons 
that will be clear later. 

The starting point for the algorithm is the same as before: computing a source vertex x using 
algorithmic The necessary LBFS ordering will be initiated from this source vertex. Just as before, 
{x}, V (G) —{x} is consistent with a linear extension and will be used as the initial ordered partition. 
Each slice will then be processed (in order) and two rounds of refinement undertaken for each. Both 
will ensure that the resulting ordered partition remains consistent with a linear extension. 

The first round of refinement targets the initial vertex of the slice in question, ensuring that 
by the end, none of its neighbours share its partition class. This has the effect of determining the 
directions on all the edges incident to that initial vertex. With that, the directions on all other 
edges active for that slice can be determined using a second round of refinement. To facilitate the 
first round of refinement, algorithm |C will be adapted so that only a single partition class (the one 
containing the initial vertex) is subdivided. 

It is critical for the running-time that during both rounds of refinement, only the edges active 
for that slice are processed. However, as will be described later, knowing the co-components for 
each slice will allow the algorithm to “cheat” a little in this regard. 

The rest of the paper is organized as follows. Section |^ outlines the various initialization steps 
that are needed, along with their correctness and running-time. Next, section |^ describes how 
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refinement can produce a linear extension, and also includes correctness and running-time. The 
paper concludes in section 


3 Initialization 

The following is needed before refinement can be applied to compute the desired linear extension: 

• A source vertex; 

• An LBFS ordering initialized from that source vertex; 

• The corresponding slices and active edges; and 

• The co-components of each slice. 

Note that a source vertex can be computed in 0{n + m) time using algorithm based on 
the partition refinement implementation of algorithm described in section |2.4[ and using that 
same implementation, algorithm is known to run in 0(n -|- m) time as well. [12] The remaining 
initialization components are addressed below in turn. 


3.1 Slices 


Explicitly computing the set of vertices in each slice would prove too costly. Instead, they will 
be implicitly computed by constructing the tree defined by the containment relationship between 
slices (see remark 2.2) as defined below: 


Definition 3.1. Let S he the set of slices defined by some LBFS ordering of the graph G. The 
corresponding slice-tree is defined such that its interior nodes correspond to the elements of S, the 
leaves correspond to the vertices in G, and the children of each slice are its maximal subslices. 


An example slice-tree is provided in figure The following are all easily verified properties of 


the slice-tree, each an extension of the containment relationship captured by remark 2.2 


Remark 3.2. Let S = 5i ,... ,Sn be the set of slices defined by some LBFS ordering and consider 
the corresponding slice-tree. 


1. If Si is an ancestor of Sj, then i < j; 

2. Si is an ancestor of Sj if and only if Sj C Si; 

3. y G Si if and only if Si is an ancestor of y. 

Recall the labels maintained by vertices during algorithm Notice that after each vertex is 
appended to the output ordering u, its label remains fixed ~ or hnal. When discussing vertices in 
the context of the slice-tree, this final label is what is meant when referring to a vertex’s label. Of 
course, slices are defined by labels, and these will be what is meant when referring to a slice’s label 
in the slice-tree context. In particular, the first slice consisting of all vertices is assigned the empty 
label. Notice the correspondence between the labels defining slices and the labels assigned to their 
initial vertex. The following are properties of the slice-tree in terms of these labels: 
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Figure 3: A slice-tree for the graph in hgure[^ in this case corresponding to the example execution 
of algorithm described in section [23| 


Lemma 3.3. Consider the LBFS ordering a and its eorresponding slice-tree. 


1. If the slice S is an ancestor of the vertex x, then S’s label is a (not necessarily proper) prefix 
of x’s label. 

2. If x and y are consecutive vertices in a, then the label of the slice being their least common 
ancestor is the longest (not necessarily proper) prefix shared by x’s label and y’s label. 


Proof. For item 1, notice that x G S' by item 3 of remark 3.2 Thus, x and S shared the same 


label at the time that S became a slice. It follows that S’s label is a (not necessarily proper) 
prefix of x’s label. For item 2, let S be the least common ancestor of x and y. Then S’s label is a 
(not necessarily proper) prefix of x’s label and y’s label by item 1; it is the longest such prefix by 
definition of S being their least common ancestor. □ 


Knowing lemma 3.3, and remembering the correspondence between vertices and slices, it is now 
easy to compute the slice-tree as in algorithm The correctness of algorithm [^follows immediately 
from the relevant definitions, lemma [3)3l and the fact that vertices are processed in order. For its 
implementation, note that identifying z is clearly a constant time operation. Obviously p can be 
computed in time on the order of the length of y’s label. And once S is identihed, the updates to 
T can be performed in constant time. 

Now consider the problem of identifying S. This can be accomplished by traversing the path 
from z to the root of the slice-tree, comparing p with the label of each slice thus encountered, 
stopping when a match is found. Of course, the last character in those two labels will be the same. 
And once those last two characters match, all the others must match as well because both are (not 
necessarily proper) prefixes of z’s label [p by definition and the other by item 1 of lemma . 
Thus, each comparison is a constant-time operation. 

The problem for the efficiency of algorithm is that several consecutive slices on the path from 
z to the root of the slice-tree might share the same label. If every label on the path was distinct, 
then identifying S would be on the order of the length of z’s label. This would mean the total cost 
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Algorithm 4: SliceTree{a) 

Input: An LBFS ordering a with initial vertex x. 

Output: The slice-tree corresponding to the set of slices defined by a. 

initialize T as the rooted tree with a single interior node having x as its first (and only) child; 
initialize the root of T by assigning it an empty label; 

foreach y ^ a — {x} (in order) do 

let z be the vertex preceding y in cr; 

let p be the longest (not necessarily proper) prefix shared by y’s label and z’s label; 
let S be the ancestor of z whose label is p; 

update T by creating a new node u and making it the last child of z; 
update T by making y the first and only child of u; 

end foreach 

return T ; 


of algorithm 1^ would be 0{n + m). The simple solution is to maintain a copy of the path from z 
to the root, excluding those nodes whose parents share their labels. Updating this path is clearly 
a constant time operation after each new vertex y is inserted. The total cost therefore remains 
0{n + m). 

3.2 Active Edges 

Active edges will drive the refinement that will later be used to compute the desired linear extension. 
Given their importance, special notation is developed for them below: 

Notation. Let 5i, ..., S'^ he the set of slices in order. For each vertex x G Si, let ai{x) be the set 
of vertices y G N{x) such that xy is active for Si. If x ^ Si, then ai{x) = 0 , by convention. 

To illustrate the notion of active edges, recall the set of slices ..., 5io defined earlier for the 
graph illustrated in figure The set of active edges for each of these slices is as follows: 


5 i: 

xb, 

by, 

bu, bz, bq, bw, br, bv, za\ 

^ 2 : 

0 ; 



53: 

yu, 

yz, 

,yq,yw,ur, qv; 

54: 

uz, 

uq. 

, UW] 

55: 

0 ; 




qw: 



57: 

0 ; 




0 ; 




12 






• ^g: 0; 

• 5io: 0. 


The slice-tree can be used to efficiently compute active edges in the manner outlined by algo¬ 
rithm below. Its correctness is based on lemma |3.4[ Notice that algorithm assigns labels to 
slices. These are not to be interpreted as the labels assigned during algorithmic Furthermore, the 
labels referenced in lemma |3.4| and corollary |3.6| afterward refer to those defined by algorithm [C 
not algorithm [C which was the custom in the previous section. 


Algorithm 5: ActiveEdges{T) 

Input: A slice-tree T corresponding to the set of slices S = Si,..., Sn defined by some 
LBFS ordering. 

Output: The slice-tree T augmented so that every slice has associated with it a graph 
induced by that slice’s active edges. 

initialize every slice Si by assigning it an empty label; 

initialize every slice Si by assigning it an empty graph referenced by G{Si); 
foreach Si (in order) do 

let y be the initial vertex for Si] 

foreach Sj,j > i, that is maximal with respect to Sj C N{y) do add y to 5j’s label 

end foreach 

foreach Sj with non-empty label do 
let Si be the parent of Sj] 

foreach pair y, z where y G Sj and z appears in Sj ’s label do 
I update G{Si) with the edge yz] 

end foreach 
end foreach 

return T ; 


Lemma 3.4. Let S be a slice and x, Si, ■ ■ ■, Sk its maximal subslices in order. If w is in the label 
for Si, then w G S — Si and Si C N{w). Conversely, for any vertex y G Si: 


1. If X G N{y), then x appears in the label of Si; 


2. If there is a z G N{y) n Sj,j < i, then z appears in the label of Si. 

Proof. It is clear from algorithm that Si C N{w). It follows that w ^ Si. Notice as well from 
algorithm IC that there is no ancestor S' of Si such that S' C N{w). Therefore, by either item 1 or 
item 2 of remark 2.4, it must happen that w G S. 

Now consider some y G Si and an ancestor S' of Si. Therefore x G S' , and so S' % N{x). Hence, 
S' is not a maximal slice such that S' C N{x). If x G N{y), then N{x) n 5 / 0, and therefore 
N{x) n 5 = Si, by item 1 of remark 2.4, So in that case, i = 1 and Si C N{x) 
will be added to Sfs label. 


It follows that X 


13 









Finally, consider some 2: G N(y) n Sj,j < i. As with x, no ancestor S' of Si can be a maximal 
slice such that S' C N{z), since 2; G S' . But Si C ^"(2;) by item 2 of remark 2.4 
will be added to Si's label. 


It follows that 2; 

□ 


Lemma 3.5. Algorithm^is correct. 

Proof. Let 5 be a slice and x, Si,..., Sk its maximal subslices in order. Consider an active edge 


xw. Then by item 1 of remark 2.4, w ^ Si. Furthermore, by item 1 of lemma 3.4, x appears in the 
label of Si. Clearly the edge xw is added to G{S) in this case. The same obviously holds for the 


active edge yz where y ^ Si and z £ Sj,i < j, although this time by item 2 of lemma 3.4 

Now consider an edge uv in G{S). According to algorithm]^ it can be assumed without loss of 
generality that u £ Si and v appears in Sfs label. Then by lemma 3.4, v £ S — Si and Si C N{v). 
It follows that uv is an active edge for 5. □ 

Consider the slices with empty labels following algorithm Let S' be a slice, x, Si,..., S^ 
its maximal subslices in order, and assume that Si has an empty label. Then by lemma [T4l the 


initial vertex x is isolated from Si as are all other vertices y £ Sj,j < i. But what about vertices 


2; G Se, In fact, no such vertex can exist, by item 3 of remark 2.4, This leads to the following 

corollary, which characterizes the relationship between active edges and the labels assigned by 
algorithm 

Corollary 3.6. Consider a slice Si and one of its maximal subslices Sj. If w is in the label for 
Sj, then for every vertex y £ Sj, w £ ai{y). Moreover, if Sj’s label is empty, then for every vertex 
y £ Sj,ai{y) = 0 . 


Proof. Follows directly from lemma 3.4 and the relevant definitions. 


□ 


In other words, the vertices of each slice without incident active edges are those in slices not 
having labels at the end of algorithm Such vertices and their containing slices are further char¬ 
acterized below. Understanding them is crucial for establishing the running-time of the transitive 
orientation algorithm developed in this paper. 


Definition 3.7. Let Si,... ,Sn be the set of slices in order. A vertex x £ Si is said to be discon¬ 
nected for Si if ai{x) = 0 and is connected for Si otherwise. 

Remark 3.8. Let S be a slice and x, Si,..., S^ its maximal subslices in order. Then G'[S'] is 
disconnected if and only if there exists a disconnected vertex for S. Moreover, when there is a 
vertex y that is disconnected for S, it happens that y £ Sk, and in this case, all vertices in Sk 
are disconnected for S. Furthermore, when S is disconnected, one of the components of G[S'] is 
{x} U U • • • U Sk-i, the others being those of G[Sk]. 

Remark |3.8| will be used below to guarantee that no disconnected vertices are processed during 
algorithmic This will be important to guarantee the linear running time of algorithm [C above. 

There are two blocks of nested for-loops in algorithm [C The case with disconnected vertices 
occurs for the second one, which will be considered below. For the first block of nested for-loops, 
notice that a simple bottom-up marking scheme can compute the necessary Sj’s for each vertex y: 
start by marking all leaves in N(y), and then recursively mark each node whose children are all 
marked. The time this takes is on the order of the size of the subtree rooted at each Sj. 
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Now, every slice contains at least one vertex: its initial vertex. Hence, if a slice contains only 
one vertex, then all of its ancestors must contain at least two vertices. Thus, the size of the subtree 
rooted at each Sj is on the order of the size of the number of its descendent leaves, or equivalently, 
the size of Sj. But the Sj^s are all disjoint, by definition, and each Sj C N{y). It follows that the 
running time of the first block of nested for-loops is 0 {n + m). 

Having finished the first block of nested for-loops, it can be assumed that the labels on each 
slice are in place. To implement the second block, it will be necessary to compute the set of vertices 
in each slice whose label is non-empty. This can be done recursively as follows. Let S' be a slice 
and X, Si, ..., Sk its maximal subslices in order. Assume that the set of vertices in each Si has been 
computed and is represented as a linked-list. For each Si whose label is non-empty, make a copy 
of the list. 

For the recursion to work, a list of all vertices in S must now be computed. Observe that 


remark 3.8 together with corollary |3.6| implies that there is at most one Si whose label is empty, 


namely Sk- If there is such an Si, take its recursively computed list of vertices and append to 
it the just computed copies of the lists from all other maximal subslices (including a copy of x). 
This produces the desired list of vertices for S. Every vertex in a slice with a label has at least 
one incident active edge, by corollary 3.6 Hence, the cost of creating the list of vertices for S is 


proportional to the total number of active edges for S. But by remark 2.3, each edge is active for 
exactly one slice. Therefore the total cost of creating these lists for all slices is 0{n + m). 

Finally, observe that once the set of labels has been computed for each slice, and the set of 
vertices has been computed for each slice with a non-empty label, the graph G{S) can clearly be 
computed in time on the order of the number of active edges for S. As above, the cost of doing so 
for all slices is 0 {n + m). 

Before moving on, it will be necessary to introduce the following results concerning active edges. 
They will be needed later to prove the correctness of the refinement that produces the desired linear 


extension. Each either rephrases or extends remark 2.4 


Remark 3.9. Let Si,..., Sn be the set of slices in order. Consider some slice Si with initial vertex 
X. If N{x) n Aj / 0, then N{x) H 5* = S'j+i = ai{x). 

Lemma 3.10. Let S be a slice and x, Si,..., Sk its maximal subslices in order. Assume that 
N{x) n 5 7^ 0. Then for every Si,i > 1, that contains a connected vertex for S, there exists a 
y G S^, i < i, that is universal to Si. 


Proof. If i ^ k, then the result follows by item 3 of remark 2.4 So assume that i = k. Suppose for 
contradiction that no such vertex y exists. Then there is no vertex in S' — that is adjacent to 


any vertex in Sk, by item 2 of remark 2.4 Therefore Si = Sk would contain a disconnected vertex, 
which contradicts remark iTSl □ 


Lemma 3.11. Let S = Si,..., Sn be the set of slices defined by some LBFS ordering. If y is the 
initial vertex for Si, then aj{y) = 0 for all j > i. 

Proof. Since y is the initial vertex of S'*, if y G Sj, then either i = j ov Sj is an ancestor of Si. 
Therefore, j < i, by item 1 of remark 3.2 So by definition, aj{y) = 0 for all j > i. □ 


3.3 Co-Components 

Section |2.6| hinted at the need to compute the co-components of the graph induced by each slice. 
To see why it is useful, consider a slice S and its maximal subslices x,Si,... ,Sk in order. Recall 
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that if iV(x) n S 7^ 0, then A^(x) n 5 = 5i, by item 1 of remark 2.4 Let C be a co-component of 


G[S'i]. Knowing the direction of an active edge incident to any vertex in C is equivalent to knowing 
the direction of the active edges incident to all vertices in C as follows: 

Lemma 3.12. Consider an LBFS ordering of a comparability graph G. Let S he one of the slices 
and X, Si, ..., Sk its maximal subslices in order. Assume that N{x) H iS / 0. Let C be a co- 
eomponent of G[Si]. Then in any transitive orientation, all edges between G and the set S — C are 
directed the same way. 

Proof. Consider some transitive orientation of G. Recall that x must be universal to G since 


N{x) n S' = by item 1 of remark 2.4 Partition G into Gout and Gin, the former being those 


vertices in G whose shared edge with x is directed toward x, the latter defined as those vertices 
in C whose shared edge with x is directed away from x. Suppose for contradiction that Gin 7 ^ 0 
and Gout / 0 - Then by virtue of G being a comparability graph, there must be a join between Gin 
and Gout, contradicting G being a co-component of G'[Si]. Since x is not adjacent to any vertex in 
V{G) — N[x], it follows that all edges between G and S2 U • • • U S^ must all be directed the same 
way as the edges between C and x. □ 

In particular, the above lemma says that once the direction on a single edge between x and G is 
known, then all the directions of the edges between x and G are known (and the directions are all 
the same). This will be useful during the first round of rehnement described in section 2.6 where 
the goal is to determine the directions on the active edges incident to x. During that round of 
refinement, the vertices in C can be moved as one block. This round of refinement therefore does 
not technically rely only on active edges. This was the “cheat” alluded to in section 2.6 


It is a simple matter to compute the co-components of each slice once the active edges for each 
vertex are known. The following result characterizes these co-components and translates directly 
into algorithmic which computes the co-components of any single slice. 

Lemma 3.13. Let S be a slice and x, Si,..., Sk its maximal subslices in order. Let Ci,..., Q be 
the co-components o/G'[S'i]. Assume that ^[5] is disconnected. Then: 

1. N{x)CS = Si; 

2. X and S 2 , ■.., Sk all belong to the same co-component of G[5], call it G; 

3. Gi C G or Gi n G = (h, for all i. 


Proof. If N{x) n 5 ^ Si, then N{x) n S' = 0, by item 1 of remark 2.4 In that case, x is a universal 


vertex in G[S], contradicting the latter being disconnected. Moreover, x must be universal to 
S2U - • • Sfc in G[S], meaning they must be in the same co-component. Now suppose for contradiction 
that Gi overlaps G. Say A = G Ci Gi and B = Gi — A. Then there is a join between A and B, 
contradicting Gi being a co-component of G[Si]. Therefore, either Cj C C or Cj n C = 0. □ 


Lemma 3.14. Algorithm^^is correct. 

Proof. If N(x) D S ^ Si, then by item 1 of remark 2.4, N{x) n S = 0. Therefore x is a universal 


vertex in GIS], and so there is a single co-component S. In this case algorithm correctly outputs 
S. So assume that N{x) (1 S = Si. Suppose that els'] is connected. Then there can be no Gi that 
is universal to C. Notice that algorithm [b correctly outputs S in this case. The case where Gl^] is 

3.13|and the relevant definitions. □ 


disconnected follows directly from lemma 
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Algorithm 6: CoComponents{S) 

Input: A slice S. 

Output: The co-components of G[S]. 

let X, Si,..., Sk he the maximal subslices of S'; 

if N(x) n S 7 ^ Si then return S 

let Cl,... ,Ci be the co-components of G[Si]; 

initialize C as the empty set; 
initialize C as {x} U S2 U • • • S^; 

foreach Ci do 

I if Ci is universal to C then C <— C U {Ci} else C •(— C U Cj 

end foreach 

c^cu{c}-, 

return C; 


All it takes now is repeated application of algorithm to compute the co-components of each 
slice. This can be done recursively as was done earlier for algorithm Let S be a slice and 
X, Si, ..., Sk its maximal subslices in order. Assume that the co-components for each Sj have been 
recursively computed and are represented as partitions of that slice using the data-structure for 
partition refinement from section [Td) 

It will be necessary to assume that for each maximal subslice, a list of its vertices has been 
recursively computed. This was also required by the implementation of algorithm and can be 
handled the same way. It will be assumed as well that algorithm precedes algorithm and hence 
the active edges for S are already known. 

Now consider algorithm operating on S in this recursive context. The maximal subslices of 
S are defined by the slice-tree. By remark 3.9 the first conditional amounts to testing if x has 
any incident active edges for S, which is a constant time operation since the active edges for each 
vertex are assumed to have been computed already. 

If the first conditional succeeds, then the list of vertices in S is returned. But this happens 
automatically as part of the maintenance required to preserve the assumption that the list of 
vertices in each maximal subslice has been recursively computed. As explained above, this can be 
handled as it is in the implementation of algorithm 

So assume that the first conditional fails in algorithm]^ The co-components of G[S'i] have been 
recursively computed, as assumed above. To initialize C, recall that the list of vertices in each max¬ 
imal subslice has been recursively computed. Copy the lists for the maximal subslices S 2 , ■ ■ ■, Sk-i, 
merge these copies into a single list, and then append the result to the list for Sk- Doing so avoids 
“touching” disconnected vertices in the case that G'[S'] is disconnected (see remark 3.8). This was 
the same approach used in the implementation of algorithm]^ (see section 3.2). For the same 
reasons as there, the total cost over all slices is 0 {n + m). 

Determining if Ci is universal to C merely requires traversing the list of active edges for S that 
are incident to each vertex in Ci- Adding Ci to C is a constant-time operation, and so too is adding 
C to C after the loop. However, adding Ci to C within the loop is on the order of the size of Ci- 
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But notice that each vertex in Ci has at least one incident active edge, namely x. Hence, adding Ci 
to C is on the order of the number of active edges for S. The cost of all this work is therefore on 
the order of the total number of active edges for S. Hence, the total cost over all slices is 0{n + m), 
by remark [2.3[ This and the preceding discussion implies that the co-components of all slices can 
be computed in time 0{n + m). 


4 Refinement 

Consider a prime comparability graph G. Based on the initialization work in section it is now 
possible to assume the following has been computed in linear-time: 

• A source vertex x G H(G); 

• an LBFS ordering with x as its initial vertex; 

• the graph induced by the active edges for each slice defined by this LBFS ordering; and 

• the set of co-components for each of these slices. 

Given the above, a linear extension for G can be computed in linear-time using partition refine¬ 
ment as described in the following subsections. 


4.1 Algorithm 

Two changes to the type of refinement defined by algorithm will be needed to produce a linear 
extension in linear-time. The hrst involves using only a subset of a vertex’s neighbours, not the 
whole neighbourhood as is done in algorithm A vertex will be pivot once for every one of their 
containing slices. The subset of neighbours processed each time will be those whose corresponding 
edge is active for that slice. One of the reasons this works relates to co-components and lemma [3.12[ 
Another is the following similar lemma: 


Lemma 4.1. Consider an LBFS ordering of a comparability graph G. Let S he one of the slices 
and X, Si,..., Sk its maximal subslices in order. If N{x)r\S 7^ 0 , then in any transitive orientation, 
all edges between St and Sj, 1 < i < j, must be directed the same way. 


Proof. Assume for contradiction that the claim does not hold for some pair Si and Sj, 1 < i < j. 
Let Si be the leftmost maximal slice for which there is a corresponding Sj. By remark 2.4 there is 
a vertex y G Si, £ < i, that is universal to Si and isolated from Sj. So by choice of 5^, if ^ > 1, then 
all of y’s incident edges with Si must be directed the same way. The same is also true if ^ = 1, by 
lemma 3.12 It follows that all edges between Si and Sj must be directed the same way: namely, 
the opposite direction of those between y and Si. □ 


Algorithmic below generalizes the pivot operation so that only a subset of a vertex’s neighbours 
are processed. Once again, the input vertex p is referred to as a pivot. Remark 4.2 is the analogue 
of remark [2T| The running-time of algorithm jC is clearly 0(|5|) based on the implementation from 
section [231 
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Algorithm 7: GeneralizedPivot{V,p, S) 

Input: An ordered partition V = Pi,..., Pk of vertices, a distinguished vertex p G Pi, and a 
set of vertices S C N{p). 

Output: A refinement of V. 

foreach Pj €'P,i^ j, such that S splits Pj do 

A^Pj- S] 

B Pj — A] 

if i < j then replace Pj in P with A, B in order else replace Pj in V with A, B in order 

end foreach 
return V] 


Remark 4.2. Suppose that A,B are partition classes (in order) that replace a partition class P 
during algorithm\^ Then A,B,S^%,p^P, and if p is in an earlier partition class than P, then 
A n S' = 0 and B C S, otherwise ACS and R n S = 0. Moreover, if P is not replaced during 
algorithm\^ then either p £ P, P C S = (), or P C S. 


The second of the required changes to algorithm involves targeting a specific partition class: 
instead of allowing all partition classes to be subdivided as in algorithm only one partition class 
specified in the input can be subdivided. However, the order of the two classes subdividing the 
original partition class must additionally be allowed to vary. This leads to algorithms and 
below. As before, the input vertex p is referred to as a pivot in both cases. Remarks 4.3 and |4.4| 
are the analogue of remark 2.1 The running-time of each algorithm is clearly 0(|S'|) based on the 


implementation from section 2.4 


Algorithm 8: PivotPush{V,p,P,S) 

luput: An ordered partition V = Pi, ... ,Pk of vertices, a distinguished vertex p £ Pi and 
partition class P = Pj, j ^ i, and a set of vertices S. 

Output: A refinement of V. 

if S splits Pj theu 
A^Pj- S; 

B Pj — A] 

end if 

if i < j then replace Pj in V with A, B in order else replace Pj in V with B, A in order 

return V] 


Remark 4.3. Suppose that A,B are partition classes (in order) that replace partition class Pj 
during algorithm^ Then either ACS and S' H R = 0 (if j > i in the algorithm) or B C S and 
A n S = 0 (if i < j in the algorithm). 
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Algorithm 9: PivotPull{V,p, P, S) 

Input: An ordered partition V = Pi,..., Pj. of vertices, a distinguished vertex p ^ Pi and 
partition class P = Pj,j ^ i, and a set of vertices S. 

Output: A refinement of V. 

if S splits Pj then 
A^Pj- 5 ; 

B Pj — A; 

end if 

if i < j then replace Pj in V with B, A in order else replace Pj in V with A, B in order 

return V] 


Remark 4.4. Suppose that A,B are partition elasses (in order) that replace partition class Pj 
during algorithm^ Then either ATS and S H B = $ (if j < i in the algorithm) or B T S and 
A n S' = 0 (if i > j in the algorithm). 


Algorithm below employs the new forms of refinement defined above to compute the desired 
linear extension. As promised in section 2.6, algorithm 10 processes each slice in order, and for 
each, there are two rounds of refinement, corresponding below to the two inner loops. 

As alluded to in section 2.6, the first round of refinement is designed to remove the neighbours 
of each initial vertex from its containing partition class. Algorithms and are used for this 
purpose, the targeted partition class being the one containing the initial vertex. During this round 
of refinement, the co-components are moved as one block for the reasons described in section 3.3 
The ordered partition maintained throughout remains consistent with a linear extension. 

Once the initial vertex is isolated from its neighbours, it can be removed from its containing 
partition class with the resulting ordered partition still being consistent with a linear extension. 
Each vertex is then processed in order and refinement takes place according to their active edges 
for that slice. With the initial vertex isolated, and by proceeding in order, this has the effect of 
refining the ordered partition in terms of the maximal subslices of the current slice. Intuitively, 
lemma [4A] ensures the result remains consistent with a linear extension. 

To illustrate the operation of algorithm 10, we return to the graph in figure and the set of 
slices 5i,..., 5io defined earlier for it in section |2.5[ It may also help to recall the corresponding 


slice tree in figure and the listing of active edges outlined in sections |3.1| and |3.2[ respectively. 

In terms of algorithm [To| we first note that V = {x}, {b, y, u, z, q, w, r, v, a} prior to any iteration 
of the outer for loop. Now consider what happens on the first iteration of the outer for loop. Note 
that each of b, y, u, z, q, w, r, v, z, and a are connected vertices for not in the same partition class 
as X, the initial vertex for Si. However, since the other endpoint of each of their active edges shares 
the same partition class as them, the inner for loop has no impact on V. Of course, neither does the 
line immediately preceding the inner for loop. The following highlights the values of V following 
each iteration of the inner for loop: 


• (After X is pivot): P = {x}, {y, u, z, q, w, r, v, a}, {&}; 

• (After 6 is pivot): P = {x},{y,u,z,q,w,r,v},{a},{b}; 
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Algorithm 10: LinearExtension{S) 

Input: The set of slices 5 = Si,..., 5,^ of a prime comparability graph G defined by some 
LBFS ordering in which the initial vertex x is a source. 

Output: A linear extension of G. 

V •<— the ordered partition {x}, V{G) — {x}; 

for Si to Sn (in order) do 

let y be the initial vertex for S*; 

let Py be the partition class in V (and throughout what follows) such that y ^ Py\ 

while there is a connected vertex z for Si that has not been pivot and such that 
z G Si — Py do 

if z G Oii{y) then 

C G- the co-component of G[ai{y)] such that z G G; 

S G- ai{z) U ai{y) - G] 

else 

S ^ ai{z); 

end if 

P G- PivotPull{V, z, Py, G) // Algorithm 
V G- PivotPush{P, z, Py, S) // Algorithm 

end while 

replace Py in P with {y}, Py — {y} in order 

foreach connected vertex z G Si (in order) do 

P G- GeneralizedPivot{P, z, ai{z)) // Algorithm 

end foreach 
end for 
return P] 
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• (After y, u, z, q, w, r, and v are pivot): V = {x}, {y, u, z, q, w, r, u}, {a}, {6}; 

• (After a is pivot): V = {x}, {z}, {y, u, q, w, r, f}, {a}, {b}. 

Now consider what happens on the second iteration of the outer for loop. Since there are no 
active edges for S 2 , there are no connected vertices for S 2 , and this iteration has no impact on V. 

On the third iteration of the outer for loop, we have the following connected vertices for S 3 : 
y,u, z,q,w,r,v, of these, z is the only one not in the same partition class as y, the initial vertex 
for S 3 . Therefore the while loop executes, and since z G 0:3(2/), the hrst branch of the if statement 
succeeds. Looking at G[a3(2/)] = G[{ri, 2;, q, rc}], we notice that C = {z, q, w} in this case; also note 
that S = {y, u, q, w}. Therefore, V = {x}, {z}, {q, w}, {y, u, r, x}, {a}, {6} after the first iteration of 
the while loop. 

Notice now that q and w are no longer in the same partition class as y. Therefore the while 

loop continues to iterate. Assume that q is selected on the next iteration. Since q G 03(7/), the 

first branch of the if statement succeeds. Once again, C = {z,q,w}, but now S = {y, x, u, z, tc}. 
Therefore, V = {x}, {z}, {q, w}, {r, x}, {y, u}, {a}, {6} after this iteration of the while loop. On the 
next iteration of the while loop, assume w is selected. Then C = {z, q, w} and S = y,u, z, q. As a 
result, there is no change, and we still have V = {x}, {z}, {q, w}, {r, x}, {y, x}, {o}, {b}. 

As a result of q and rx’s refinement on the previous iterations of the while loop, we no longer 
have r and x sharing the same partition class. Therefore the while loop continues to operate. 
Assuming r is chosen next, the second branch of the if-statement succeeds because r ^ 0:3(2/)- 
Hence, 0 = 0 and S = a3(r) = {x}. Therefore, V = {x}, {z}, {q, w}, {r, x}, {y}, {x}, {a}, { 6 } after 
this iteration. Assuming now that x is chosen on the next iteration, we have 0 = 0 (because 
the second branch of the if statement succeeds) and S = {g'}, meaning there is no change and 
r = {x},{z},{q,w},{r,v},{y},{u},{a},{b}. 

It is now no longer possible to further refine y’s partition class. Thus, even though x is a 
connected vertex not in t/’s partition class, when it is chosen, which will have to be on the next 
iteration of the while loop, it will have no effect on V. After this, all connected vertices for S 3 not 
in 2/’s partition class will have been selected and no further iterations of the while loop will execute. 

Moving on, observe that the line immediately prior to the inner for loop also has no impact on 
V due to y being in a singleton class. We now consider what happens to V on each iteration of the 
inner for loop. 

Recall that the connected vertices for 53 (in order) are: y,u, z,q,w,r,v. The following demon¬ 
strates the state of V after each of these vertices is selected as pivot on subsequent iterations of the 
inner for loop: 

• (After y is pivot): V = {x}, {z}, {q, w}, {r, x}, {y}, {x}, {a}, {b}; 

• (After X is pivot): V = {x}, {z}, {q, w}, {r}, {x}, {y}, {x}, {a}, {b}; 

• (After z, q, w and r are pivot): V = {x}, {z}, {q, w}, {r}, {x}, {2/}, {x}, {a}, { 6 }; 

• (After X is pivot): V = {x}, {z}, {g}, {x;}, {r}, {x}, {2/}, {x}, {a}, {b}. 

Notice that all partition classes are now singleton sets. As algorithm clearly never unions any 
two partition classes, we can stop our tracing of the algorithm as V will not be further modified. 
Observe that V defines a linear extension of G corresponding to the directions on the edges in 
figure 
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4.2 Implementation and Running-Time 

4.2.1 Preprocessing 

First consider the portion of algorithm [T^ outside the outer for loop. Recall from section that 
we can assume we have access to the underlying LBFS ordering defining the slices with x as the 
initial vertex. We can therefore assume access to x in 0{n) time. Given x, the initial partition 
{x},F(G) — {x} can be formed in 0 {n) time using the partition refinement data-structure from 
section [2^ Before moving on to the outer for loop itself, we will need to perform some preprocessing 
to ensure its efficient execution. 

First, we need to compute the initial vertex for each slice. But recall the correspondence 
between slices and vertices. The initial vertex for Si is therefore the vertex in the LBFS 
ordering. Therefore the initial vertex for each slice can be computed in 0{n) time. Given this, 
access to Py can be obtained in constant time throughout the outer for loop using a simple array 
data structure. 

Next, we will need to compute an ordered list of connected vertices for each slice, where the 
ordering is determined by the LBFS. This will specifically be needed for the efficient execution 
of the inner for loop. To help us, we’ll use the graph induced by the active edges for each slice, 
which we computed earlier during initialization (see section]^; we’ll assume this is represented as 
an adjacency list. Notice that it suffices to order the vertices in this adjacency list according to the 
underlying LBFS ordering. 

To do this, we’ll need to first make a copy of the LBFS ordering, to which as we noted above, 
we can assume access. Then, we’ll need to associate a linked list (initialized as empty) with each 
vertex in that copy. After doing this, we’ll process each vertex in the graph induced by the active 
edges for Si, appending a reference to itself in the list associated with it in the LBFS copy. We’ll 
then do the same for each of S'2 ,..., 5'n in order. Afterwards, the list associated with each vertex in 
the LBFS copy clearly references its position in the adjacency list of every graph induced by active 
edges (as long as it appears there). In other words, the list corresponds to those slices where the 
vertex is connected. 

The next step is to go through the copy of the LBFS in order, traversing the list associated 
with each vertex in order, following each reference, then removing that vertex from the graph’s 
adjacency list only to then append it to that same adjacency list. After doing this, it should be 
clear that the list of vertices for every graph induced by the active edges for each slice has been 
ordered according to the underlying LBFS ordering. 

The cost of doing this is clearly proportional to the sum (over all vertices) of the number of 
slices each vertex belongs to for which it is connected. By definition, each vertex is connected for a 
slice if it has an incident active edge for that slice. But each edge incident to a vertex is active for 
at most one slice, by remark [2.3[ Hence, the total cost to order the vertices in the graph induced 
by the active edges of each slice is 0 {n + m). 

4.2.2 The Outer For Loop 

Now consider each iteration of the outer for-loop. As we noted above, each initial vertex y and 
its partition class Py can be located in const ant-time. Notice that access to Py is only needed a 
constant number of times on each iteration. Hence, the total cost of this access over the entire 
algorithm is 0 {n). 
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To evaluate the condition in the while loop, we will maintain a list of connected vertices for Si 
that have not been pivot, and we will partition this list into two halves: those that are in Py and 
those that are not. Given such a list, the condition in the while loop can clearly be evaluated in 


constant time using the data structure of section 2.4 


Such a list will first need to be constructed prior to any iteration of the while loop. This can be 
performed by scanning the list of vertices in the graph induced by the active edges for Si, to which 
we have access as noted earlier. At each vertex, it is a constant time operation to verify it is in Py 
assuming a pointer to Py and the data structure for partition classes from section |2.4[ The total 
cost of initializing this list for each Si over the entire algorithm is therefore 0(n + m), by applying 


remark 2.3 in the same fashion we did above. Maintaining this list of vertices (that has not been 


pivot) is clearly proportional to the time to maintain V itself, and hence we turn to that now. 

Recall from section that the co-components of each Si can be assumed to be computed. 
Therefore, if the first branch of the if-statement in the while loop succeeds, the vertices in C can be 
identified and removed from Py in \C\ time during the call to algorithm]^ A key observation for 
the running time is that this only needs to be done once. That is, once a vertex from C has been 
pivot all vertices in C will have been removed from Py already. Therefore, on any iteration of the 
outer for loop, the vertices in C only need to be processed once. Of course, each of these vertices is 
adjacent to the initial vertex for the slice being processed, and each of the incident edges is active 
for that slice. So once again, by remark 2.3, the total cost of algorithmis 0{n + m). 

Consider now the total cost of identifying and removing S from Py over the course of the 
algorithm. Of course, S can be computed and removed from Py by traversing ai{z) and ai{y), 
to which we have constant time access via the graph induced by the active edges for Si, as noted 
earlier. But that might potentially lead to traversing ai{y) multiple times during the while loop, 
which would not be consistent with linear time overall. 

The solution is to process ai{y) and y as one “block”. What we mean by this is that the first 
time ai{y) is processed, the portion in Py should be made children of the vertex y in the partition 
refinement data structure of section 2.4 Then, for the duration of the while loop, ignore the 


requirement to include ai{y) as part of S. This works because the only time when ai{y) must be 
added to S is when the first branch of the conditional in the while loop succeeds, which happens 
when y G ai{z): if the vertices of ai{y) H Py are children of Py, then they will move with y as part 
of ai{z) as they should. 

It may happen however that some portion of ai{z) n Py needs to be included in S when the 
second branch of the conditional succeeds. In this case, the vertices in question can be removed 
from “under” y as though y were a nested partition class within Py. This adds at most a constant 


amount of work to each step in the partition refinement implementation of section 2.4, At the end 
of the while loop, one more traversal of ai{y) can restore Py to its original condition. 

Overall, then, the identification and removal of S from Py results in ai{y) being traversed at 
most twice. The remaining work involves traversing ai{z) on each iteration of the while loop. This 


work is consistent with linear time overall 0(n -|- m) due to remark 2.3 as above. 

Since the line immediately preceding the inner for loop is clearly a constant time operation, it 
only remains to consider the inner for loop itself. Using the partition refinement implementation of 


section 2.4 the cost of each iteration is clearly proportional to ai{z). Each vertex is only pivot once, 
so overall, the cost of the inner for loop is clearly consistent with linear time, meaning 0 {n + m). 
From this and our work above, we can conclude the following: 


Lemma 4.5. Algorithm 10 can he implemented to run in 0{n + m) time. 


24 










4.3 Correctness 


Having shown that algorithm p!0| can be implemented to run in linear-time, it only remains to prove 
its correctness. To this end, notice that on each iteration of the outer for-loop, the initial vertex 
of the slice under consideration is isolated in its own partition class. Thus, the ordered partition 
at the conclusion of the outer for-loop induces an ordering of the vertices. It therefore suffices to 
prove that the ordered partition maintained throughout algorithm 10 is consistent with a linear 
extension. In order to accomplish this, the corresponding invariant needs to be strengthened as 
below. 

Invariant 4.6. After the iteration of the outer for-loop: 


1. The endpoints of any edge that is active for the sliee Si reside in different partition elasses in 
V, and if i > 0, eaeh endpoint of such an edge is universal to the partition class containing 
the other endpoint; 

2. If i < n, then there is at least one active edge with respeet to Si+i whose endpoints reside in 
different partition classes; 

3. The ordered partition V is consistent with a linear extension. 


The invariant provides some intuition as to the operation and correctness of the algorithm. 
Condition 1 says that at the end of each iteration of the outer for-loop, no further refinement is 
possible using active edges with respect to already considered slices. Moreover, condition 3 says 
that all refinement up until that point has been consistent with a linear extension. And finally, 
condition 2 says that on the next iteration, there will be at least one “seed” edge to begin the 
process of refining using the active edges with respect to the next slice under consideration. 

Starting from this seed edge, refinement takes place in a manner and order that ensures that all 
active edges for that slice will be used to refine the ordered partition maintained by the algorithm, 
and in such a way that is consistent with a linear extension. By the end of the algorithm, every edge 
will have been used to refine the ordered partition. All such refinement will have been consistent 
with a linear extension. As remarked above, only singleton sets will remain, thereby inducing a 
linear extension. 

The rest of this section constitutes one large induction proof of invariant |4.6[ Each part - base 
case, induction hypothesis, induction step - will be stated and proved separately. All statements 
should be interpreted in the context of algorithm [Tol that is, G, Si,... ,Sn, x, y, Py, C, S, etc. 


have the meaning assigned to them in the algorithm. 


Base Case. Invariant f.6 holds before any iteration of the outer for-loop. 


Proof. The initial ordered partition {x}, V{G) — {x} is consistent with a linear extension by virtue 
of X being a source vertex. Since the input graph G is prime, it must be connected, and therefore 
X must have at least one incident active edge with respect to the slice Si, the other endpoint of 
which is obviously in a different partition class than is x. The remaining condition of invariant |4.6| 
holds vacuously. □ 


Induction Hypothesis. Consider the iteration of the outer for-loop. 
after all prior iterations. 


Then invariant f.G holds 
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Induction Step This will be comprised of three separate lemmas, one for each condition of 
invariant 14.61 


Lemma 4.7. Assuming the induction hypothesis, condition 1 of invariant 4-6 holds after the i^^ 
iteration of the outer for-loop. 


Proof. To prove the first part of condition 1 of invariant 4.6 we will instead prove the following 
claim: immediately prior to the iteration of the inner for-loop, all active edges with respect to 
Si that are incident to the first j + 1 vertices in Si are such that their endpoints reside in different 
partition classes in V. 

This cearly holds prior to any iteration of the inner for loop because the initial vertex for Si 
resides in its own partition class prior to any iteration of the inner for-loop. So consider the 
iteration of the inner for-loop and assume the claim holds after all prior iterations. 

Let z be the j vertex in Si and choose P G V he such that z G P. Suppose for contradiction 
that there exists a vertex v G ai{z) such that v G P as well. Notice that v cannot be amongst the 
first j -|- 1 vertices in Si. Given the base case, it will be assumed that z is not the initial vertex for 


Si. Therefore, by condition 3 of remark 2.4, there is a vertex u that is earlier than z in Si such 
that u is adjacent to z but not adjacent to v, and in this context, z G ai{u) while v ^ ai{u). 

Suppose that u is the vertex in Si. Then i < j. So by the induction hypothesis, immediately 
prior to the i— 1'^* iteration of the inner for-loop, u and z resided in different partition classes in V. 
However, v and z must reside in the same partition class at this point as partition classes are clearly 
never merged during the inner for-loop. So by remark 4.2, immediately prior to the iteration of 
the inner for-loop, v and z must reside in different partition classes in V. The desired contradiction 
follows by once again noting that no union ever occurs between a subset of one partition class and 
a subset of another partition class during the inner for-loop. 


This proves the claim and implies condition 1 of invariant 4.6 Note that the second part of that 


condition follows as a corollary to the claim by applying remark 4.2 to every non-vacuous iteration 
of the inner for loop. □ 


Lemma 4.8. Assuming the induction hypothesis, eondition 2 of invariant 4-6 holds after the i^^ 
iteration of the outer for-loop. 

Proof. If there is no active edge for Sj+i then the lemma follows vacuously. So assume that at least 
one such edge exists. Note that both its endpoints are connected. Now assume for contradiction 
that after the iteration of the outer for-loop, the endpoints of all such active edges reside in 
the same partition classes in V. Let C be the set of connected vertices for Si+i. Note that C" is 
non-trivial by assumption. Furthermore, there is a partition class P G V such that C C P. Notice 


as well that C is a non-trivial module in G[S'j+i], by remark 3.8 But since G is prime, it is not 
a non-trivial module for G, hence, there is a vertex q ^ S'j+i and vertices u,v G G' such that q is 
adjacent to precisely one of them, say u. Let Sj be the least common ancestor of q and u in the 
slice tree. Note that j < i -|- 1, by item 1 of remark 3.2, Note as well that uq is an active edge 


for Sj, by choice of Sj. So by condition 1 of invariant 4.6 and our induction hypothesis, u and q 


reside in different partition classes at the end of the iteration of the outer for loop. Note that 
at no point during algorithm |10| is there any union between some subset of one partition class and 
a subset of a different partition class. Therefore u and v are in the same partition class after the 
jth if; 0 ];-a,tion of the outer for loop. So by condition 1 of invariant 4.6 and our induction hypothesis, 
q must be adjacent to v, a contradiction. □ 
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Lemma 4.9. Assuming the induction hypothesis, eondition 3 of invariant 4-6 holds after the i^^ 
iteration of the outer for-loop. 

Proof. Consider the i + 1^* iteration of the outer for loop. We will prove that V is consistent with 
a linear extension after this iteration by proving the following three claims; 

1. "P is consistent with a linear extension after every iteration of the while loop; 

2. N(y) (1 Py = ^ after the while loop has completed; 

3. V is consistent with a linear extension after every iteration of the inner for loop. 


Claim 1: Notice that V is consistent with a linear extension prior to any iteration of the while 
loop by our induction hypothesis. So assume V is consistent with a linear extension after some 
number of iterations of the while loop. Note that on any iteration of the while loop, V only changes 
by the splitting of Py as a result of calls to algorithm and algorithm 

After the call to algorithm P is replaced in V with Py — C,C H Py (not necessarily in that 
order). Observe that if V is no longer consistent with a linear extension, it is because of an edge 


with one endpoint in — C and the other in C H Py. However, by lemma 3.12, these edges must 
share the same direction as those between z and C, and therefore V remains consistent with a 
linear extension after the call to algorithm 

For the call to algorithm]^ Py is replaced with Py — S and P^nS (not necessarily in that order). 
Observe that 2 is universal to Py O S. We argue that all edges between Py — S and PyCiS must be 
directed the same as those between Py — (iS and z. Suppose for contradiction there is an edge that 
is not, between say u G Py — S and v G Py (1 S. If z is not adjacent to u, then the contradiction 
follows immediately. Therefore u and z must be adjacent. 

We now have two cases. First, assume that z G ai+i{y). Then u ^ because if it were, 

then either u G aj+i(y) or u G aj+i(z), neither of which hold by definition of S in this case. So let 
Sj be the least common ancestor of u and v in the slice-tree. Notice that j < i + \ and the edge uv 


is active for Sj. Also notice that since algorithm 10 clearly never unions a subset of one partition 


class with a subset of another partition class, and since u,v G Py prior to the call to algorithm 
it follows that u and v were in the same partition class at the conclusion of the iteration of the 


outer for loop. But this contradicts condition 1 of invariant 4.6 and our induction hypothesis. 

So now assume that z ^ aj+i(?/). Then S = ai+i(z), meaning z and v are in different maximal 
sub slices of Sj+i. Furthermore, u ^ Q:i+i(z), meaning either u ^ S'j+i or there is a maximal subslice 
of Si+i, call it S, such that u,z G S. A similar argument as in the previous case rules out the 


former possibility. However, the latter case is ruled out by lemma 2.5 


Claim 2: Let M be the set defined as the union of {y}, A^(y)nPynSj+i, and all maximal subsets 
of that are subsets of Py. We will show that M is a module for G[5i+i]. Consider any vertex 
q G Si+i — Py. We have two cases, depending on whether q is adjacent to y. First consider the 
case where q is adjacent to y. Then by definition, q G ai+i{y), and so q is connected for Sj+i. As 


q ^ Py, we know that q was pivot on some iteration of the while loop. So by remark 4.3 and the 


definition of S in algorithm 10 in this case, q is universal to Py and hence universal to M . 


Now consider the case where q is not adjacent to y. So in this case q ^ aj+i(y). Let S' be the 
maximal subslice of Sj+i such that q G S'. Note that there is no vertex u G S' H M by definition of 
M and since q ^ Py. Also note that if q is disconnected, then so are all vertices in S' , and in this 
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case, all vertices in S' are isolated from M, by definition. So assume that q is connected. Then it 


was pivot on some iteration of the while loop, since q ^ Py. So by remark 4.3 and the definition of 
S in algorithir|10| in this case, we have that q is isolated from Py and hence M. We will now show 
that all vertices in (S' H Py) — M are similarly isolated from M to conclude that M is a module for 
G[Si+. 

Consider some vertex r G (S' Ci Py) — M. Assume for contradiction that r is adjacent to some 
vertex u G M. Notice that u cannot appear in an earlier maximal subslice than S' due to item 2 
So u appears in a later maximal subslice than S', call it S". Then by item 3 
there is a vertex v ^ S' U S" that is universal to S' and isolated from S". So 


2.4 


2.4 


of remark 
of remark 

q,r € ai+i(v) but u ^ Q;j+i(u), and v is connected for Sj+i. And recall that q is isolated from Py. 
Therefore v ^ Py. Hence, q was pivot on some iteration of the while loop. Since S" appears later 
than S', and q ^ ai+i(y), we also have u ^ ai+i(y). Thus, after the iteration on which v was pivot. 


we cannot have u,r ^ Py as we do, by remark 4.3 


lemma 4.8 


From above, it follows that M is a module for G[S'j+i]. Applying a similar argument as in 
we can conclude that M is a module for G. However, since G is prime, this module 
must be trivial. But we know y G M, by definition, so it must be that N(y) n Py n S'j+i = 0. In 
other words, no neighbour of y in appears in Py. We conclude the claim by citing item 1 of 
invariant 14.61 and lemma 13.111 


Claim 3: By claim 2 above, V is consistent with a linear extension prior to the first iteration 
of the inner for loop. So assume it is consistent with a linear extension after some number of 
iterations of the inner for loop. Consider the pivot ^ as defined by the next iteration of the inner 
for loop. Suppose that during this iteration, 2 splits a partition class, replacing it with A,B, in 
order according to algorithmNote that z ^ Au B. Assume that z appears in a partition class 
earlier than A in P (the other case is symmetric). It suffices to show that all edges between A and 
B are directed toward B. 

Let uv be such an edge, and say u £ A and b £ B. It follows that u ^ ai+i(z) and v £ ai+i(z). 
If u ^ A'(z), then the edge uv must clearly be directed toward v as it is here. So assume that 


u £ N(z). Then by condition 1 of invariant 4.6 along with the induction hypothesis, it follows that 


the edge uz is not active for any Sj,j < i + i. Therefore u and 2 must be in the same maximal 
subslice for By lemma 2.5, it suffices to consider the case where u,z £ N(y). Here we have 

u,z £ ai-\.i(y) as well. 

Observe that v ^ y because on the prior iteration of the inner for loop, both u and v resided in 
the same partition class, and clearly y resides in its own singleton partition class throughout the 
inner for loop. Since v £ ai+i(z), it follows that v ^ N(y) in this case, and therefore v ^ ai+i(y). 
It follows that on the first iteration of the inner for loop, y would have split the partition class 
containing u and v, contradicting our assumption regarding A and B. It follows that u,v ^ P^(y)- 
Therefore V is consistent with a linear extension after each iteration of the inner for loop. □ 


5 Conclusion 

This paper describes a linear-time algorithm to compute a linear extension of a prime compara¬ 
bility graph. It does so by extending the elegant partition refinement framework of [TTj. Some 
initialization is required as part of its implementation, but this work only requires tree and list 
traversals. This represents a significant simplihcation over the only other linear-time algorithm for 
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the problem due to m- It was achieved by using LBFS to accomplish the same effect as modular 
decomposition in the earlier algorithm. 

The main problem with that earlier algorithm was its reliance on the notoriously difficult linear¬ 
time modular decomposition algorithm from that same paper. A simpler, linear-time alternative 
has since been developed in [23] and could likely be used in its place, but it would not accomplish 
the same degree of simplification as observed here with LBFS. Nevertheless, a reduction from 
modular decomposition is still required to compute a linear extension for non-prime comparability 
graphs, and for that purpose, the simpler, linear-time modular decomposition algorithm of |23j 
is the preferred choice. It can be viewed as the companion to this paper. Many of its insights 
contributed to the development of the algorithm presented here. 

It is important to draw the distinction between the transitive orientation problem and the 
comparability graph recognition problem. As already stated, the transitive orientation problem asks 
you to orient the edges of a comparability graph. In that case, you are given a comparability graph. 
In contrast, the recognition problem gives you an arbitrary graph and asks you to determine if it is 
a comparability graph. One solution would be to apply the transitive orientation algorithm from 
this paper and then verify if the resulting orientation is actually transitive. Rather surprisingly, 
this verification is the harder problem, at least in terms of running-time. The book [22| has a good 
discussion on this. As snmmarized there, the problem has been shown to be equivalent to finding a 
triangle in a graph, and in both cases, the fastest algorithm requires matrix multiplication, whose 
complexity is currently |6|. Some kind of hardness result in terms of matrix multiplication 

or any more efficient algorithm would be a significant result. 

Fortunately, a transitive orientation of a comparability graph is frequently all that is required 
for many applications M- Other applications involving co-comparability graphs and permutation 
graphs require a transitive orientation of the complement graph. One of these applications is 
permutation graph recognition. For efficiency reasons, it is preferable if an orientation of the 
complement can be computed without explicitly computing that complement. This is possible by 
implicitly representing such an orientation by a linear extension. With minor modihcations, the 
previous linear-time transitive orientation algorithm of m can do so in linear-time. The algorithm 
developed in this paper is capable of the same. Among other implications, this would result in a 
simpler linear-time permutation graph recognition algorithm than the one presented in |15j . which 
is the only other to date. 

Two changes are required to transform the algorithm developed in this paper so that it computes 
a transitive orientation of the complement graph. First, as in [T2|, the order of the subdivided classes 
in algorithm would need to be reversed. Doing so has the effect of computing a source vertex 
for the complement while also maintaining consistency with a linear extension of the complement 
throughout rehnement. Moreover, [12] showed that the same change results in an LBFS ordering 
of the complement when applied to the partition refinement implementation of LBFS. The second 
modification concerns the special role played by the co-components of slices. This would now 
have to be played by the components of slices. But to avoid “touching” disconnected vertices and 
thereby maintain linear-time, rather than co-components being “pulled” out of Py by their other 
halves during the first round of refinement, components will have to be “pushed out” from within 
Py to join their other halves. The proof of correctness and running-time would be similar to the 
one in this paper. 
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